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Abstract 

We propose a phenomenological equation for the vortex line density in rotating Bose- 
Einstein condensates as a function of the angular speed. This equation provides a simple 
description of the gross features of the increase in vortex number from the appearance of 
the first vortex to the theoretical rigid-body result for high vortex density, and allows one 
to compare with analogous situations in superfluid helium, after the suitable changes in 
the relevant parameters are made. 

1 Introduction 

The study of quantized vorticity in superfluid helium 4 has been a relevant topic in superfluidity 
for several decades H]. In the last decade, much interest has been focused also on the study 
of quantized vortices in rotating Bose-Einstein condensates of alkali-atomic gases [2]-flU], for 
several geometries of the conflning potential. Here we focus our attention on the vortex density 
L per unit transversal area in rotating BEC in terms of the angular velocity 17, between the 
lower critical velocity for the appearance of the flrst vortex to fast rotations, where vortex 
discreteness may be neglected, and the system is expected to be characterized by the rigid 
body result L = 2Q/k, with k = h/m being the vorticity quantum, h Planck's constant and 
m the mass of the particles. 

To compare the behavior of vortices in both systems, helium 4 and BEC, it is interesting 
to show macroscopic analogies and differences arising from the respective microscopic features 
of these systems, as for instance the low density, the long coherence length, the size of the 
vortex radius, and the compressibility of atomic BEC as compared with the high density, the 
extremely thin vortices, and the incompressibility of liquid helium 11. Furthermore, due to very 
low temperatures, all particles of BEC participate in the superfluid component, whereas in 
liquid helium only a part of the particles are in the superfluid component and the other ones 
in the normal viscous component. For the mentioned comparison, it may be useful to explore 
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phenomenological frameworks common to both systems in which the macroscopic behavior 
may be put in a common ground. 

The aim of this paper is to adapt an evolution equation for the vortex hne density proposed 
for rotating superfluid hehum [111 [T2] to the vortex density in BEC superfluids, by taking 
into account the physical differences in the respective situations. This is done in Section 
2, whereas Section 3 is devoted to the study of the solutions of the equation and of their 
stability, from which the vortex density is obtained in terms of the angular velocity and the 
corresponding vortex number is calculated. In Section 4 we discuss a physical interpretation 
of the macroscopical equation. In Section 5 we discuss a special situation where the number 
of vortices seems to have a maximum instead of being a monotonically increasing function of 
the angular velocity. 



2 Evolution equation for vortex density 



Recently, two of us proposed a phenomenological equation for the evolution of the vortex line 
density L per unit volume in superfluid helium 4 under rotation and counterflow [TT], [T^]. In 
the case of pure rotation in a cylindrical container, such an equation takes the form 



d 



^3/2 



(2.1) 



where d is the diameter of the container, and and (3i dimensionless phenomenological coef- 
ficients, whose values are obtained by fitting experimental data. 

Here, we aim to explore the application of a suitably modified form of this equation to 
the description of the vortex density in rotating BEC in terms of the rotation frequency. The 
consequences of equation ()2.ip have been studied in detail for superfluid ^He [Ej, where the 
situation without vortices L = becomes unstable for values of higher than a critical value 
^1 = (/33/2/32)^'^/'^^) and the vortex density L increases for increasing Q beyond the critical 
value f^i. For high values of $7, it follows from (|2.ip that the steady-state solution L exhibits 
the well-known behavior L = 2VI/k [1]. Thus, equation (j2.ip describes the increase from one 
single vortex to many vortices, and it is more suitable than taking directly the assumption 
L = 2Vt/K. Note that for parallel vortices of the same height, L is simply the number of 
vortices per unit transversal area, i.e. L = N/A, with A the total transversal area of the 
system. Therefore, in this situation the vortex line density per unit volume coincides with the 
vortex line density per unit transversal area, to which we will pay attention here for rotating 
BEC. 

In contrast with superfluid helium, Bose-Einstein condensates are not confined in material 
cylindrical containers, but in optical or magnetic potential traps, to which one may add a 
magnetic rotation in analogy with rotating container experiments for superfluid helium. It 
is often assumed that the trap has the axially symmetric form V{r) = ^mO.ir'^ with r the 
transversal radius and fit the frequency trap, which is a constant characterizing the confining 
potential, or an ellipsoidal form V{x,y,z) = m{uj^x'^ + w^y^ -|- J\z^)l2^ where lOx^ and 
ujz are the oscillator frequencies characterizing the potential of the trap in the three spatial 
directions, and that the X and Y axes rotate with a frequency fixed by the frequency of 
variation of a magnetic field along the X and Y axes [6]. 

Comparing with the rigid cylindrical containers in usual experiments on rotating ^He su- 
perfluids, in Bose-Einstein condensates one has not a deflnite value for the diameter d of the 
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container. Here we will take instead of it the characteristic size of the trap, which controls the 
spatial extent of the one-particle ground state of the harmonic potential and which is called 

1 /2 

the oscillator length, aosc ~ {k/uj±) ' . It is also useful to recall, for future applications, that 
in axially symmetric traps the stationary number particle density in terms of the radius in the 
Thomas- Fermi approximation is [71 [TO] 



n(r) 



h 



-L0± 



1/2 



P2 



(2.2) 



with N the total number of particles of the condensate, a the s-wave scattering length appearing 
in the Gross-Pitaevskii equation [2]-[TU], Iz the average extent of the condensate on the z 
direction — taken as the direction of the axis of rotation — , lo±_ the average transversal 
frequency of the trap, given by uj\_ = (w^. -|- u;y)/2, or simply by Vlt in symmetric traps, and 
Rtf is the Thomas-Fermi radius of the condensate, given by [Tj llOj 



Rtf = 2 



Na 



1/4 



1/2 



To apply equation (|2.ip to rotating BEG, we identify d as Coscj and (|2.ip becomes 



dL „ ^2 / 

— = -pKL + a2VKW± 



a2 



L, 



(2.3) 



(2.4) 



where = i}/uj± is the dimensionless angular velocity. A further relevant difference with 
the usual rotating cylinder used in superfluid helium may be the anisotropy of the trap. In 
this case the coefficients of (12. 4|) become a function of the eccentricity e, defined as e = {cOy — 
io'^)/{uj'^+(jjy). We do not pretend that equation ()2.4p yields a full explanation for the dynamics 
of the BEG, which would require to include into the description other collective modes, but 
we aim to identify some of the most salient phenomenological analogies and differences with 
superfluid helium, which should, in the future, to be understood from a microscopic basis of 
the Gross-Pitaevskii equation. 



3 Stationary solutions of the evolution equation 

Now, we study the stationary solutions of equation (|2.4p and their corresponding domain of 
stability. The stationary solutions are 



L = 0, 
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(3.5) 
(3.6) 



The first solution (L = 0) corresponds to the absence of vortices, and the second one 
describes an increase in the vortex density for increasing Q. Regarding the stability of solution 
L = (no vortices), and according to equation ()2.4p . we consider the evolution equation for 
the perturbation 6L around L = 0, which is 



5L 
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(3.7) 
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From this equation we can estabhsh that the solution L = is stable if 



Pin-P3Vn + a4 > 0. 

In superfluid helium [12], a satisfactory choice of the parameters is to assume /3| 
which case the critical value of $7 is 



2/3i 



a4 



(3.8) 

4a4/3i, in 

(3.9) 



and L = is stable for Q < Jlc. For symmetric traps, the critical value in BEC is = 1/V2, 
corresponding to half the frequency of the quadripole mode [5]-[9] which is a collective surface 
mode with angular momentum I = 2; when resonantly excited, quantized vortices formed on 
the surface come into the condensate. In eccentric traps, the value of flc depends on the 
eccentricity e, and it is close to Q{e) = — e (some authors [T] give this value whereas other 

ones [6] indicate Q{e) ^ ^-0.91e). In any case, the critical value of O for the formation of the 
first vortices is lower in eccentric traps than in symmetric traps, which reflects the fact that 
the finite energy barrier that vortices must overcome to move into the condensate is lowered 
if the condensate is elongated. 

For Q > r^c, the behavior of L in terms of Q is described by ()3.6p . which tends to L = 
2i}L0±/K for high values of fi. 
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Figure 1: Vortex density L in a BEC in terms of the dimensionless angular velocity for two 
values of the eccentricity e = (continuous line) and e = 0.037 (dashed line). The graphics 
are plotted using the value w_L/27r = 219 Hz for the the average transversal frequency of the 
trap from Ref. [3]. 

Under the previous assumption = 4a4/3i, equation ()2.4p may be rewritten as 

dL n t1 , 

— = -PkL + a2^JKUJl_ 



(3.10) 
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whose stable solution for 17 > f^c is simply 
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(3.11) 



The value of the combination of coefficients appearing in the prefactor of (j3.1ip is dictated by 
the fact that for 17 S> tends to the rigid body result L = 2^}uj±/k. In Fig. [T] we plot L in 

terms of $7 for a given uj± and for two values of eccentricity of the potential, e = (continuous 
line) and e = 0.037 (dashed line) [1]. As expected, graphics of L in figure [T] confirm the 
fact that when the eccentricity increases the first vortex in BEC appears for a smaller value 
of rotation. Expression ()3.1ip generalizes the well-known expression L = 20, /k, and, from a 
practical perspective, it is more useful than it, because in BEC, in contrast to helium 4, the 
angular velocity cannot be increased indefinitely, because of the effects of the centrifugal force. 
Note, of course, that these lines are not completely realistic, and only discontinuous values of 
them, corresponding to integer numbers of vortices, are actually meaningful. 
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Figure 2: The number of vortex Ny in a BEC in terms of the dimensionless angular velocity 
for two values of the eccentricity e = (continuous line) and e = 0.037 (dashed line). The 
graphics are plotted using the values a;_L/27r = 219 Hz, A'^ = 10^, a = 8 fim and Iz = 49 nm 
from Ref . [3j . The horizontal lines indicates the ranges of to which the number of vortices 
corresponds: the first range corresponds to 1 vortex while the second one corresponds to 2/4 
vortices. 

In many analyses, one measures the number Ny of the vortices rather than their areal 
density. To go from L to Ny one must multiply L times the transversal area of the condensate, 
taking into account that the value of its radius depends on the rotation rate, because of the 
centrifugal force [LQl- One has 

Nv = LttK = — "2 = =2"' v-^-l^) 

with Rtf the Thomas- Fermi radius ()2.3p . Note, then, that for a given value of the angular 
frequency the total number of vortices increases when the number of particles increases, pro- 
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portionally to (Na/lz) 



thus, it depends not directly on N, but rather on the number of 



particles per unit transverse length. A second aspect is that cannot increase indefinitely, 
because at = uj± the centrifugal force is so high that the condensate is flattened out and 
it may disperse, unless a suitable quartic potential is added to the harmonic trap. However, 
since both L and R increase with Q, it follows that as far as the condensate exists, the total 
number of vortices given in ()3.12p increase as function of the angular velocity, as shown in 
Fig. [21 where Ny is plotted for two values of the eccentricity e = (continuous line) and 
e = 0.037 (dashed line). There a direct comparison between the graphics of our model and 
the experimental data of Ref. [1] is carried out. Note, incidentally, that the relation between 
L and Ny is very different in BEC than in confined superfluid helium, where the number of 
vortex is 



with R the radius of the cylinder, k the vorticity quantum — which is different with respect 
to that of BEC — and J7c the critical rotation velocity for the appearing of the first vortex 
in rotating helium II. The difference between two relations ( (|3.12|) and (|3.13|) 1 is caused by 
the fact that in the latter the radius of the sample always coincides with the radius of the 
container and it does not depend on the angular velocity, whereas in BEC the change in the 
radius with is very relevant. 

To compare the behavior of vortices in rotating Bose-Einstein condensates and in rotating 
superfluid helium 4, it is also worth to stress that the characteristic values of the angular 
velocities and of the vortex densities are rather different in both cases; in typical experiments 
with helium the angular frequency is less than 10 rad/s, whereas in BEC may be of the 
order of 300 rad/s; and the vortex line densities are of the order of 10^ vortex lines/cm^ for 
superfluid helium and of 10'' vortex lines/cm^ in BEC. However, we should not compare the 
actual values of 17 but characteristic dimensionless rotation velocity. In BEC, it is defined as 
r2 = r2/tJ_L, whereas for superfluid helium it is given by (7 = QcP/k; for ^He, the value of the 
quantum of vorticity k is k ~ 9.97 x 10~^cm^s~^; thus, the value of the dimensionless velocity 
corresponding to 10 rad/s for helium in a container of diameter 1cm, would be of the order of 
10^, much higher than those explored in Bose-Einstein condensates. This difference can also 
be related to the fact that the quanta of vorticity of helium II and BEC are very different, as 
the quantum of vorticity in rubidium BEC is about 20 times lower than in helium II because 
of the lower atomic mass of helium with respect to that of rubidium. 

The equation (j2.4p has been applied to superfluid helium in rotating cylinders. It would be 
interesting to have information for rotating containers with ellipsoidal cross section, in order 
to explore the influence of the eccentricity on the vortex line density, but we are not aware 
of experimental results in this problem. A parallel possibility would be to compare rotating 
containers with rectangular cross section |13j . 

4 A physical interpretation of the evolution equation 

To have a more explicit understanding of equation ()2.4p or (I3.10p we recall that in rotating 
Bose-Einstein condensates and Helium II, the vortices are produced at the surface of the con- 
densate and penetrate into it pulled by the rotating drive. The repulsive interaction amongst 
them tends to push them apart because vortices rotating in the same direction experience 




(3.13) 
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an effective repulsive interaction [T, "5]. This competition between centripetal driving and re- 
pulsive force yields eventually to a regular vortex lattice. This idea suggests us to identify 
the terms appearing in equation (12. 4p in terms of these features. In this way, we propose to 
interpret (12. 4p as 

^ = [Vin - Vout]L^^\ (4.14) 

with Vin, Vout being the respective ingoing and outgoing vortex drift velocities, proportional 
to the respective driving forces. The exponent 3/2 for L in (I4.14|) is chosen on dimensional 
grounds, since L has dimensions of (lenght)"^. There is not an univocal way to decide to which 
term, either Vin or Vout-, must be assigned each term in (12. 4p or in or (I3.10p . Here we propose to 
relate the terms linear in L^/^ to the inflowing flux, because the corresponding velocity would 
not depend on the vortex density, but only on the rotation. Following this tentative criterion, 
the respective values of Vin and Vout would be 

2 

(4.15) 



The outgoing velocity depends on the vortex line density: it depends on the inverse of 
the average vortex separation (which is given by L^^/^). Thus, (|4.15p could be considered as 
the first terms in a kind of virial expansion of a repulsive force between vortices in terms of 
powers of L^/^. In vortex arrays in rotating containers, the free energy has a term which may 
be interpreted as a repulsion force. The ingoing velocity is related to the angular velocity: 
it is positive at values of i7 higher than Vtc- This could probably be related to an energy 
barrier which the vortices must overcome in order to move into the condensate. Then, the 
interpretation of a dynamical equation for L for rotating superfluids — BEC or liquid helium 
— is rather different than for counterflow systems, where it relies on the dynamics of vortex 
breaking and reconnections, as proposed Schwarz [13], whose microscopic view is not useful 
for parallel vortex lines. Furthermore, in Schwarz microscopical model for the dynamics of 
L under counterflow in superfluid helium, all the volume contributes to the production term 
of L, because of the elongation of vortex loops, whereas in rotating systems the vortices are 
produced on the surface. A relevant difference between BEC and superfluid helium may be the 
extent of vortex interaction, because of the considerable width of vortices in BEC as compared 
with the atomic width of vortex lines in superfluid helium. This effect should be reflected in 
the value of coefficient (3. 



Vout = I3kL^'^ + 13 




5 An anomalous non-monotonic density variation 

The usually expected behavior of L in terms of 17 is a monotonical increase as mentioned in 
Section 3. However, in an experiment by Hodby et al. [6] a different behavior was observed, 
which we briefly report here for the sake of completeness, although it seems rather exceptional, 
as it has not been reproduced, to our knowledge. In [6] Hodby et al. studied experimentally 
the evolution of the number of vortices for a condensate of 2 • 10"^ atoms of ^"^Rb for an oblate 
geometry {ijo±_ < uJz) within a rotating magnetic elliptical trap. They observed that when the 
eccentricity is adiabatically ramped from e = to a given final value e, there were both a lower 
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r2i(e) and an upper (e) values of rotation rates for which vortices were nucleated; these 
values depend on the eccentricity, and were given respectively by [6] 



Oi(e) « 0.71 -0.91e, 



and 



2 



2Qo 



2/3 



(5.17) 



(5.18) 



Hodby et al. [6] measured the number of vortices as a function of O for two values of the 

within the interval f^i < 17 < 



, by identifying its coefficients in terms 



eccentricity, and found that this has a maximum {L. 
^2- This information may also be synthesized into 

of r2l(e), rJ2(e), ^max{e) and Lmax{^)- 

To describe Hodby's results one should assume /3| > Aa^Pi instead of the equality of both 
terms assumed in (|3.9|) . In this case, there will be an instability range for L = for 0, between 
r^l^^^ and ^l^^"^ given by 



2/3i 



4q4/3i 



2/3i 



a 



1/2 ^ P3_ 



+ 



2/3i 



(5.19) 



(5.20) 



from which we deduce immediately that f^^/Pi = + \/Q2 and 04/(^1 
Now, we rewrite the non zero stationary solution (|3.6|) as 



"2V'^-L 

2^ 



1 ± 



1 + 



4/3i 



/2 



(5.21) 



with 02 = Ci2/(3 and /3i = We have taken 03/02 = ■\/r22 in order to ensure that the 

solution of L will vanish at J7 = 1^2. This condition implies a restriction on the coefficients in 
([23D, namely 03/02 = P3/P1 - (o402)/(/5i03). _ _ 

The branch is not considered here because inside the interval [r2i,J72] it is negative; 
whereas -L_ vanishes at f^i and $72 and it has a maximum for a value of the frequency V^l^rnax ? 
given by the condition 



02 



VPi VO2 - 

The corresponding 



0,2 ~ V^r 



2 • 



(5.22) 



(5.23) 



Thus, the coefficients in (12. 4p have been identified in terms of Oi{e), 1^2 (e)j ^max{^) and 
Lmaxi^), from which their values could be found, except for /3, related to the dynamics of the 
vortex line density. The coefficients Oj, Pi depend on the eccentricity e, but not on uj± neither 
on the mass of the particles, which appears in k. 
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To compare with the experimental values of Hodby et al. [6] one must go from the values 
of L (number of vortices per unit area) to the average number of vortices in the condensate, N , 
which are related through L = N/A, with A the transversal area of the condensate. This will 
be estimated as ^ = {'7r/4)d? = {t: /A){k/u}±). In fact, it must be recalled that the transverse 
radius depends on the angular velocity O as commented in ()3.12p . In the trap used by Hodby 
et al. [6], uj± ~ 1247rs~"'^, and for ^''Rb, the vorticity quantum k = h/m is 4.54 x lO^^cm^s"""^. 
Thus, we take for the average area of the condensate A ~ 9.15 x 10~^cm^. In Fig. 3 the curves 
of our model and the experimental results by Hodby et al. are shown [6]. The experimental 
data are disperse but the general qualitative features are well exhibited. We have taken for Oi, 

^max and Lmax their experimental values — in fact, the experimental value for 0.2 found 
by Hodby et al. does not coincide exactly with the theoretical value coming from (j5.18p . 

Hodby et al. give the value of ^max for two values of e, namely, e = 0.041, for which 
^max ~ 0.74, and e = 0.084 for which Q. 

max ~ 0.785. In view of the values reported by these 
authors, we suggest for ^rnax{^) the ansatz 

^max{e) « ^£±M) = ^ + 102(6). (5.24) 

In fact, the values of 1^2 (e) observed in Fig. 1 of [6] were r22(e = 0.041) ~ 0.77 it 0.02 and 
^2(e = 0.084) 0.86 ± 0.02, in such a way that (I5.24p yields indeed a reasonable estimate for 
^max{^) in these cases. 

6 Conclusions 

We adapted a previous dynamical equation proposed for the vortex density in superfluid helium 
to BEC condensates. The adaptation is not obvious, because it requires to consider the form 
of the trap potential, the radius of the BEC cloud in terms of the number of particles and the 
angular velocity, and information on the critical value of Q related to the resonant excitation 
of the quadrupolar mode. 

An interesting aspect of our proposal is the simplicity of the expression (jS.lip for the vortex 
density in term of the angular rotation. The lines shown in the figure are indicative, as in fact 
the number of vortices only take discontinuous values, whereas the actual values are related 
to the discrete number of vortices. Relation (jS.lip shows the considerable difference between 
the vortex density starting from L = at = ■^^i. and tending to L = ^ for an asymptotic 
limit of i7 ^ L0±_ which is never reached, as the centrifugal force limits O to values less or 
equal than a;_|_, with uj±_ the transversal frequencies of the confining trap. Thus, (j3.1ip albeit 
phenomenological, may be useful for simple estimates of the density, and also of the number 
of vortices, as given by (|3.12p . 

We also included in Section 5 a seemingly anomalous situation where the number of vortices 
is not an increasing function of Vt, but is reported to have a maximum. To our knowledge, 
this observation has not been replicated by other groups, and it could be related to some 
metastable situation. In any case, we have reported it here to have a better appreciation of 
the role of some restrictions on the values of the numerical coefficients, as that reported in 
()3.9p . whose breakdown could lead to a situation as that reported in Section 5 instead to the 
much more generic results of Section 3. 

A physical interpretation of the meaning of the several terms in ()2.4p is still open, and is 
related to the entrance into the condensate of vortices formed in the walls, and their subse- 
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quent mutual interactions inside the condensate. 
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Figure 3: Representation of the number of vortices N in terms of the dimensionless rotational 
velocity according to equation (|5.21|) . The points correspond to the experimental data of 
Hodby et al. [6]. The continuous lines correspond to our model, taking for VLi, 0,2, ^max and 
Lmax the experimental values. The upper line corresponds to the eccentricity ei = 0.084 and 
the lower one to eccentricity 62 = 0.041. 
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